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Abstract 

The folding transition of single, long semiflexible polymers was studied with spe- 
cial emphasis on the chain length effect using Monte Carlo simulations. While a 
relatively short chain (10-25 Kuhn segments) undergoes a large discrete transition 
between swollen coil and compact toroid conformations, a long chain (50 Kuhn 
segments) exhibits an intrachain segregated state between the disordered coil and 
ordered toroid. 



1 Introduction 



The conformational transition of a polymer chain from a statistical coil to 
a condensed state has been actively studied over the past several decades. 
According to the classical scenario, a polymer coil gradually shrinks into a 
spherical globule with a decrease in the solvent quality or an increase in the 
pair-wise attraction between monomers [1]. However, once the effect of the 
chain stiffness is taken into account, the situation changes, even qualitatively. 
For a single semiflexible polymer, the collapse transition becomes discontinu- 
ous and shows the characteristics of a disordered disperse state and an ordered 
condensed state p]l2l[3lll] . This transition has been demonstrated by the single 
molecule observation of giant DNA molecules 
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Under usual aqueous conditions, DNA possesses a rigid helical structure, which 
gives rise to a large persistence length Ip ~ 50 nm compared to its molecular 
thickness c? ~ 2 nm. Therefore, giant DNAs with contour length L S> 1 /xm is 
characterized as a typical semifiexible polymer chain. There have been many 
numerical and theoretical studies on the folding transition of a semifiexible 
chain [Tll3|l^|l5|l6|l7|8ll9l[T0|ll|12|13yi4|15|16 These studies show a variety of 
collapsed morphologies, such as a toroid, a rod, and a spherical globule, along 
with strong finite chain length effect. It has been shown that a toroid is the 
most typical structure as the compact state, whereas some other morphologies 
are also found, such as a rod, a sphere, and a composite between a toroid and 
rod. Although observations by computer simulation have been shown to be 
very powerful, the range of chain length studied so far has been rather lim- 
ited; i.e., the length is too short, which is at most on the order of 10 Kuhn 
segments. On the other hand, genomic DNA molecules have a contour length 
on the order of mm~cm, indicating that the number of Kuhn segments is on 
the order of 10^-10^. Even for phage DNAs, the number is on the order of 10^. 
Thus, studies on long semifiexible chain would have biological significance, in 
addition to physical interest p!8]. On the other hand, theoretical studies usu- 
ally postulate some specific morphologies and analyze their relative stabilities 
by a grand state or mean-filed type approximation. Overall, our understanding 
of the collapse transition of long semifiexible polymers is far from complete, 
and it would be desirable to describe a generic scenario that could be expected 
for longer chains. 

In this Letter, we report the nontrivial scenario of the collapse transition found 
in long semifiexible polymers based on a systematic Monte Carlo simulation. 



2 Methods 



To investigate the collapse transition of long semifiexible polymers, we carried 
out off-lattice Monte Carlo simulations with a Metropolis algorithm in three- 
dimensional space, in which the polymer is represented as N beads (diameter 
a) connected by springs. The potential energy of the system is represented by 
the following three terms; 
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where Fi is the coordinate of the ith monomer, 9i is the angle between adja- 
cent bond vectors and T is the thermal energy with the Boltzmann constant 
kB = 1. The length and energy are measured in units of a and e, respectively. 
We set the spring constant /cbond = 40 and impose a cut-off, i.e., Ub^ads = oo, 
if I Fi+i — Ti |> l.lScr or | ri+i — Fi |< 0.85(7 to avoid unrealistic bond stretch- 
ing. For the bending elasticity, we choose kg = 60, which corresponds to a 
constant persistence length Ip ~ 10a independent of T. The excluded-volume 
and attractive effects are included in the Lennard- Jones potential Uij. We set 
N — 200, 500, and 1000, which correspond to a DNA chain of about 3, 7.5, 
and 15 kilo base pairs, respectively. 

For the trial motion, we adopted a slithering-snake motion, in which an end 
monomer is removed and attached to the other end of the chain, as well as the 
usual random-hopping method. The introduction of a slithering-snake motion 
makes sampling in the configurational space very efficient and enables us to 
obtain statistics even for long chains. This is particularly useful in transition 
regions, which allows for transitions among different (meta) stable states sev- 
eral times in a single long simulation on the order of 10^ ~ 10^ Monte Carlo 
steps. Each calculation was started from a randomly generated swollen coil 
conformation, and the statistical quantities were sampled every 5000 Monte 
Carlo steps after equilibration. 

We introduce two "order-parameters" concerning the statistical quantities 
here. One is the quantity of the folded fraction of monomers, and is defined 

as 



where is an indicator of the pair contact: pi^j = 1 if the number of 
monomers, which satisfies | Fi — rj |< rd, is more than 1, and otherwise 
pij = 0. Therefore, p{i) represents a binary local density (p = 1 if the i-th 
monomer has some surrounding neighbors and p = otherwise). The other is 
following local orientational order 77, 



77 = ^(3(C0S^ Oij)iocal - 1) (5) 



where the bracket {)iocai indicates averaging over the bond pairs \ i — j \> X, 
the spatial distance between which is less than rc2- In the following discussion 
we set Tci = 2.5(7, rc2 — 1-5(7 and A = 3. 
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3 Results 



Figure 1 (a) shows the average of the normahzed gyration radius {Rg)/N^ 
as a function of the inverse temperature. The bracket (•) denotes the ensem- 
ble average. A swollen coil collapses into condensed states with a decrease in 
temperature, and the transition is discontinuous, as shown by the bimodal 
distribution of Rg around the transition point. To clarify this point, we divide 
a total ensemble into subensembles with different states (here the coil and con- 
densed states) and calculate the average quantity in each subensemble. Figure 
1 (b) shows the fluctuations of the gyration radius {{Rg) — {Rg)'^)/N. A large 
peak is observed in the transition region, which reflects the co-existence of the 
coil and condensed states. These data indicate that with an increase in the 
chain length, the transition temperature becomes higher and the discontinuity 
increases. 

Figure 1 (c) shows the fluctuation of enthalpy {E"^) — {E)"^, where E = 
Ubond + Ubend + ULj, to chcck the chaugc in properties upon a folding transition 
beside the parameters that are related to the gyration radius. For short chains 
{N — 500 and 200), the temperature dependence showed only a single peak 
at e/T ~ 0.47 and e/T ~ 0.56, respectively, as shown in Fig. 1 (c). The single 
peak corresponds to a large discrete transition between the coil and compact 
states, where the latter state is characterized by the morphology of a compact 
toroid, which is consistent with previous studies. In contrast, for a long chain 
{N = 1000), a pair of peaks are noted at e/T ~ 0.41 and e/T ~ 0.44, indicat- 
ing that the folding transition proceeds via two steps. 

Three-dimensional snapshots are shown in Fig. 2 (a) at e/T = 0.46 (after the 
second peak) and in Fig. 2 (b) at e/T = 0.43 (between the first and second 
peaks). While Fig. 2 (a) is a regular toroidal structure, the long scmiflcxible 
chain also exhibits a qualitatively different condensed state: a core-shell struc- 
ture (Fig. 2 (b)), which consists of a condensed core and a coiled shell around 
it. The overall shape of the core is spherical, and is topologically different from 
a toroid. In contrast to a toroid, there is no angle to describe the features of 
a core-shell. With a decrease in temperature, the first transition of the coiled 
state to a core-shell is in the first peak, e/T ~ 0.41, and the second transition 
of a core-shell to a toroid is in the second peak e/T ~ 0.44. When e/T ~ 0.43, 
the initially prepared toroid structure spontaneously transforms into the core- 
shell structure, which suggests that the core-shell is a globally stable structure 
in this temperature range. 

Figure 3 (A) shows the distribution of rj at three different temperatures for 
N — 1000. The degree of local orientational order clearly shows a difference 
between the toroid and core-shell. In the core-shell (e/T = 0.43) and toroid 

(e/T = 0.46) states, single peaks are seen near t] = 0.2 and t] = 0.4, respec- 
tively. On the other hand, two or three peaks are observed at e/T = 0.44, 
which indicates the co-existence of structures with different orientational or- 
ders. 
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Fig. 1. (a): Normalized gyration radius (Rg) versus inverse temperature e/T. In 
the transition region, where the distribution becomes bimodal, each subensemble 
is treated separately (see the text). The lines are simply visual guides, (b): Fluc- 



tuations of gyration radius 



{Rg)^)/N versus inverse temperature e/T. (c): 



Fluctuations of enthalpy {E ) — (E) versus inverse temperature e/T. 

Figure 3 (B) shows the peaks of orientational order rjpeak versus inverse tem- 
perature. For = 1000 (Fig. 3 (a)), the folding transition is characterized 
by a two-step mechanism; the first step corresponds to the collapse from a 
swollen coil to a core-shell, and the second involves the transition from a core- 
shell to an ordered toroidal structure with a decrease in temperature. In the 
state when the core-shell and toroid co-exist {e/T ~ 0.445), a close inspection 
reveals that there are three states: the core-shell, the toroid, and another with 
a corresponding peak in local orientational order at around r^peafc — 0.3. Based 
on snapshots, the structures would be a core-shell of different type with not a 
spherical but a toroidal core, or a disk, which is similar to a toroid but there is 
no hole in the center. The data with N = 500 (Fig. 3 (b)) indicate the possible 
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Fig. 2. Typical snapshots of the collapsed structure for a chain with N = 1000. (A) 
a toroid at e/T = 0.46 and (B) a core-shell at e/T = 0.43. The bar represents the 
persistence length; Ip = 10a. 




occurrence of structures with rjpeak — 0.3 even at this chain length within a 
narrow temperature interval (perhaps as metastable states). 
Figure 4 (a) shows the bending energy per unit {Ubend)/N. While the toroidal 
structure is favorable in terms of bending, the core-shell has a larger bend- 
ing energy than a toroid and a coil. Note also the small peak in the case of 
= 500, which again indicates the appearance of core-shell or disk structures 
in this region. Figure 4 (b) shows the average fraction of folded segments 
(P). (P) of the core-shell is clearly lower than that of the toroidal structure 
((P) ~ 1), which indicates the partially folded state, and thus the cohesive 
energy of the core-shell is lower than that of the toroidal structure. These re- 
sults suggest that the stability of the core-shell has an entropic origin because 
the core-shell has a greater enthalpy than the toroidal structure. 



6 



(A) 



cr 






0.12 


£= 




o 




ribut 


0.08 






T3 


0.04 


Q) 




\- 










core-shell — 
co-existence 
toroid — 



-0.4 




0.8 1.2 



(B) 



(a) 



(b) 



0.6 - 

^ 0.4 - 
o 

I 0.2 

Q. 

- 



(0) 



0.2 0.3 0.4 0.5 0.6 0.7 

8/T 



0.6 

^ 0.4 - 
o 

I 0.2 - 

- 



0.2 0.3 0.4 0.5 0.6 0.7 

6/T 




0.2 0.3 0.4 0.5 0.6 0.7 
e/T 



Fig. 3. (A) Distribution of the local orientational order rj for a chain with N = 1000 
at various temperatures e/T. Solid, dashed and dotted curves correspond to 
e/T = 0.43 (core-shell), e/T = 0.44 (co-existence) and e/T = 0.46 (toroid), re- 
spectively. (B) Peak of orientational order rypeafc versus inverse temperature e/T. In 
coiled state, there are no peak in rj, so the averages of orientational order (rj) are 
plotted, (a): For N = 1000, there are four states; coil (rj) ~ 0, core-shell r/peafc — 0.2, 
core-shell (with toroidal core) or disk fjpeak — 0.3, and toroid fjpeak 

~ 0.4 ~ 0.6. 

(b): For N = 500, there are three states; coil, core-shell (with toroidal core) at 
Vpeak — 0.2 ~ 0.3, and toroid. (c): For = 200, there are two states; coil and 
toroid. 



4 Discussion 



The stability of the toroidal structure has been discussed in detail by several 
authors [TT|14fl5fl6fl7j . These theories claim that the toroid becomes fat with 
a decrease in solvent quality and/or upon lengthening of the chain, which leads 
to shrinkage of the toroid hole. As a result, the toroid is assumed to transform 
into a disk or sphere as low temperature structures. A rod-like structure has 
also been reported to be stable in some parameter regions [5|T0] . All of these 
studies considered the stability of a condensed state at sufficiently low temper- 
atures, i.e., well below the collapse transition point, where all monomers show 



7 



N = 1000- 
N = 200 ■■ 



N = 500 - - ■ 



(a) 



V 




0.2 0.3 0.4 0.5 0.6 0.7 

e/T 



(b) 



1 . 
0.8 
0.6 
0.4 
0.2 




0.2 0.3 0.4 0.5 0.6 0.7 

e/T 



Fig. 4. (a): Average bending energy (f/feend) versus inverse temperature e/T. (b): 
Fraction of folded segments {P) versus inverse temperature e/T. 

the condensed state (-P = 1). While various morphologies have been suggested 
due to the finiteness of the chain length, they are implicitly assumed to be in 
the same ordered "phase". 

In contrast to the current understanding, in this study we demonstrated that 
a single long semiflexible polymer can assume distinct condensed "phases" 
characterized by different orientational orders. Such a chain can also assume 
a partially collapsed state with P < 1. These results show that the structure 
of a single long semiflexible polymer can be tuned upon the folding transition. 
These features may be relevant to the biological function of DNA molecules 
through the control of higher-order structures. 

T4 DNA (166 kilo base pairs, about 500 Kuhn segments) and A- DNA (48 kilo 
base pairs, about 150 Kuhn segments) are examples of giant DNA that are 
often used in experiments, and it is expected that they are long enough to ex- 
hibit a core-shell structure. There are some experimental indications (though 
not conclusive) regarding the appearance of the core-shell structure as a col- 
lapsed state of these giant DNA molecules near the transition region. First, 
fluorescent microscopic observations often reveal "blurred" images of collapsed 
DNAs around the transition region, which is signified by the fluorescent inten- 
sity distribution as well as relatively slow Brownian motion. Second, distinct 
elastic responses in an experiment on the mechanical unfolding of A- DNA 
have suggested the presence of different condensed states [I9] . While the stick- 
release pattern observed in the deep condensed region most likely reflects the 
ordered toroid structure, another pattern with a long plateau in the force- 
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distance curve which is characteristic of the shallow condensed region may 
be attributed to a more loose and disordered structure. Furthermore, as more 
direct evidence, electron microscopic observations have shown that a single T4 
DNA molecule assumes a partially folded core-shell structure and/or a loosely 
packed folded structure under certain conditions [20] . 

The large discrete nature of the folding transition of giant DNA has been 
demonstrated in single DNA observation through the use of fluorescence mi- 
croscopy. Although such observations have clarified this discrete nature, the 
detailed microscopic structure of the product of the folding transition has only 
been visualized using methodologies such as electron microscopy and atomic 
force microscopy. Unfortunately, these methodologies only provide information 
on species adsorbed on a solid surface. However, some experiments indicate 
the presence of a loose, less ordered collapsed state compared to an ordered 
toroid in the shallow condensed condition, and it may be reasonable to inter- 
pret these features in terms of the core-shell structure found in the present 
study. 



5 Conclusion 

We studied the nature of the folding transition of a semifiexible polymer chain 
with changes in the contour length through Monte Carlo simulation. For rel- 
atively short chains (10 or 25 Kuhn segments), this transition is all-or-none 
type and produces a toroidal structure. In contrast, for long polymer chains 
(50 Kuhn segments), this transition is a multiple-step process, in which a 
swollen coil is first collapsed into a partially folded core-shell structure, and is 
followed by the subsequent transitions to the completely folded ordered struc- 
ture upon further quenching. Our findings, and in particular the appearance 
of a partially folded core-shell as a stable structure, may be unexpected based 
on the current understanding in this field. A statistical analysis indicated that 
the core-shell is an entropically stabilized state, and its appearance should 
be rather ubiquitous in the folding of long semifiexible chains. The fact that 
long semifiexible polymers possess such structural variability without a drastic 
change in the spatial size may have important consequences in the context of 
DNA function in vivo. Further theoretical and experimental studies are needed 
to fully elucidate the folding transition of long semifiexible polymers. 
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